Phase Structures of Magnetic Impurity Models with Two-Body Hybridization 
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The most general model with a magnetic impurity coupled to hybridizing and screening channels 
of a conduction band is considered. The partition function of the system is asymptotically equivalent 
to that of the multi-component kink plasma with a weak external field. The scaling properties of the 
models for finite U are sketched by using the Anderson- Yuval-Hamann-Cardy poor man's scaling 
theory. We point out that it is proper to include a two-body hybridization in order to obtain 
correct renormalization flows. The phase structures are studied graphically for the general model 
and various reduced models. A Fermi-non-Fermi liquid phase transition is found for all the models. 
We also show all possible phases with different finite temperature behaviors though they have the 
same Fermi liquid fixed point at low temperature. We also discuss the fixed point behaviors in the 
mixed valence state regime. 



I. INTRODUCTION 



The normal-state properties of the Cu-O-based ma- 
terials had been understood phenomenologically on the 
basis of the existence of the non-Fermi-liquid(NFL)-like 
groundstatcs Jl[] Q. In one dimension (Td), it is well- 
known that the Luttinger liquid states appear univer- 
sally in a wide class of interacting systems [p| . The NFL 
behaviors are shown for arbitrarily weak interactions 
in 1-d. However, for the dimensions higher than one, 
the perturbative renormalization group analysis does not 
support the NFL behaviors Q. This promotes the 
non-perturbative approaches to be developed. In 1-d, 
bosonization is a powerful non-perturbative method ||. 
For higher dimensions, several ways are applied. People 
try to generalize the 1-d bosonization procedure to higher 
dimensional one. Haldane ||, Houghton and Marston |t]] 
provided a kind of bosonization procedures in higher di- 
mensions. More recently, many authors use this kind 
of procedure to discuss whether the NFL behaviors ap- 
pear in higher dimensions ||. The other subject that 
many investigations focus on is the confinement of the 
Luttinger liquid behavior among the weakly interacting 
1-d Tomonaga-Luttinger chains |J . There still are some 
arguments for this subject |10|. The way related to this 
paper is to study such models that the properties of the 
models can be obtained either exactly or with the well- 
defined limit of validity. For example, various single im- 
purity models are good objects to be studied. Research 
for them can be thought of as a starting point to do the 
corresponding lattice models related to the physics of the 
high temperature superconduction materials. 

It has been known that at low energies which are much 
less than the Kondo temperature, the groundstates of 
simple magnetic impurity models ]i"l| ] Jl2| are Fermi liq- 
uid (FL) [[13) juj. It has recently been argued that the 
generalization of those models is closely related to the 
physics of the Cu-O-based materials @ @ Q @. It 
is important that a phase transition between FL and NFL 
states has been found in investigations of single impu- 
rity models because the NFL behavior can explain many 
anomalies of the normal state of high temperature super- 



conductors § §. 

The Kondo model is one of typical impurity models. 
Through the bosonization or the partition function ana- 
log, the Kondo problem can be mapped to a one-channel 
resonant-level model and has been shown that the ground 
state for the energy well below the Kondo temperature is 
of the FL behavior jl9). The multi-channel generalized 
model considers all of conduction channels interact with 
the impurity orbital but only one channel hybridizes with 
the impurity orbital. Following the terminology used by 
the authors of |H| jl6) we call later the 0-channel or 
the hybridizing channel and others the screening chan- 
nels. For the finite range interactions a FL-NFL phase 
transition has been shown by multiplicative renormal- 
ization grou p |2C[| i n [lq| and by the poor man's scal- 
ing theory plj in |]I7|. Physically, as pointed out by 
Zhang, Yu and Su |j!7||, the phase transition is caused 
by a competition between two different kinds of metallic 
behaviors, the NFL behavior of the multi-channel X-ray 
edge problem in the screening channels and the FL 
behavior of the one-channel Kondo problem in the hy- 
bridizing channel. A more general impurity model is the 
Anderson model jll| [i2fl . Perakis et al |l5| studied the 
multi-channel Anderson model and also found that there 
are NFL phases by using Wilson's numerical renormal- 
ization group method |13[| . A similar case also occurs for 
a generalized Hubbard model in infinite dimensions [ p"8| . 

In a previous Brief Report by one of us, co- 
operating with Li and d'Amberumenil, one showed that it 
is necessary to introduce a hybridization of two hybridiz- 
ing electrons with opposite spin hopping onto impurity 
orbital in order to keep the renormalization group's tran- 
sitivity in the multi-channel models. Here, we would like 
to give our full investigation to this issue. We will gen- 
eralize the approach used in the particle-holc(PH) sym- 
metric model to more general magnetic impurity mod- 
els in the spirit of the Anderson- Yuval-Hamann-Cardy 
(AYHC) poor man's scaling theory |2l[] ^4j, which has 
been generalized to a form with symmetry breaking [ p5[ . 
We also find that in the most general model for finite U 
a fully renormalizable model must include a hybridiza- 
tion of two hybridizing electrons with opposite spin hop- 
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ping onto impurity orbital. Because we try to describe 
the low energy behaviors under a low electron density 
in the conduction band, we apply the bosonization pro- 
cedure for impurity models ||26f . The Hilbert space of 
the system can be divided into four sub-spaces character- 
ized by the impurity states. The 'dipole' operators be- 
tween sub-spaces are of their transitivity. This transitiv- 
ity leads to an asymptotic equivalence, via bosonization, 
between the general model and a multi-component kink 
plasma model. This has been shown by Si and Kotliar 
in the infinite U case [ p~8|j25| ] . On the other hand, ac- 
cording to the transitivity, the 'charges' of the plasma 
are shown to obey the Jacobi identities, which implies 
that the 'charges' form an adjoint representation of some 
Lie algebra. The Lie algebra is shown to be solvable. 
Under the PH symmetry, the dimensions of the Lie alge- 
bra are reduced, which is closely related to the changes 
of the phase structures. The phase structures of various 
models can be recognized from the corresponding renor- 
malization group equations which are derived by using 
the AYHC poor man's scaling theory and its symme- 
try breaking generalization. The phase diagrams for the 
models may be drawn systematically by Cardy's graphi- 
cal method p3 . For each model, we find that the phase 
in which all fugacities of the plasma are irrelevant is the 
NFL phase. As long as one of the fugacities varies from 
an irrelevant to a relevant one, the FL-NFL phase transi- 
tion occurs. For the irrelevant phase, the model at fixed 
point is reduced to the X-ray edge problem, whose fixed 
point behaviors are controlled by the NFL. A compe- 
tition between those two different type metallic states 
causes the FL-NFL phase transition. Cardy's graphical 
method also shows that it is possible that there are phase 
transitions at finite temperature between other phases al- 
though they have the same low temperature Fermi liquid 
fixed point. However, the renormalization group equa- 
tions may not be valid quantitatively in the strongly cou- 
pling regions. We also study fixed point behaviors in the 
mixed valence regime. One sees that there are four kinds 
of the different fixed points. One is in the strong cou- 
pling limit corresponding to the FL one and another is 
in the weak coupling limit which is the NFL one. We 
calculate various single-particle and two-particle corre- 
lation functions in the weak coupling limit. Other two 
are the intermediate coupling regimes in which there are 
local NFL behaviors although the infrared orthogonal- 
ity catastrophe of the X-ray edge type does not appear. 
Comparing to the results of the infinite U model |25|] , we 
have one more intermediate mixed valence regime which 
is related to the two-body hybridization irrelevant. 

This paper is organized as follows. In section II, we 
define the general model under the consideration and ex- 
plain various reduced models. In section III, the par- 
tition function of the theory is expressed by the histo- 
ries of the impurity states and is mapped into that of 
a multi-component kink plasma. We also briefly review 
the bosonization procedure for the problem. In section 
IV, we study the phase structures of the theory. In sec- 



tion V, the fixed point behaviors of the mixed valence 
states are discussed. The last section gives our conclu- 
sion. The appendices talk about the Lie algebra related 
to the 'charges' of the kinks and give some comments on 
the derivation of the renormalization group equations. 



II. SINGLE IMPURITY MODELS WITH 
TWO-BODY HYBRIDIZATION 



We begin with the derivation of the general model. As- 
sume a single impurity atom is located at the origin and 
has a localized electron orbital (f>d(r) with the local im- 
purity level tda- Other electrons are in the conduction 
band described by the wave function ^k(r) and energy 
£k- A general state function is the summation over all 
possible states, 

V^(r) = 0k(r)X CT c kCT + </>d(r)X da d a , (2.1) 

ker a 

where X a and Xd a are the spin wave functions and a 
labels spin up \ and down J,. More generally, we consider 
the spin wave function of impurity, Xd a as an operator 
Xda = x a {\ — nd) + x a nd, which leads to the different 
wave functions for the empty, single occupation and dou- 
ble occupation of the impurity orbital. Under this con- 
sideration, the Fano- Anderson model is generalized as 



H G fa = ^ e k c L c k<r + tdad\d a 



(2.2) 



k.tr 



—?f y^[(*kq(l - n d ) + t^ndY^da + h.c], 
V L 



k.tr 



where tda is the impurity level which is taken as the spin 
independence, e<j, without losing the generality. Assume 
only the conduction electrons in the 0-channel hybridize 
with the impurity orbital and the hybridizations are k- 
independent in the model. Th e ch emical potential is set 
to zero but the Hamiltonian ( |2.2| ) is still particle-hole- 
asymmetry. 

To discuss the interesting physical phenomenon, the 
electron-electron interaction must be introduced. We 
consider only the finite range interaction. Besides the 
Hubbard-type interaction in the Anderson model we also 
keep consideration to other two types of interactions, the 
spin-exchange interaction and a two-body hybridization 
which describes the two hybridizing electrons with op- 
posite spins hopping onto the impurity orb ital simulta- 
neously. By using the wave function (2.1), the model 
Hamiltonian reads 



H = ^ e kc\ al c ka i + e d n d + U n d ^n d i 



(2.3) 



k>0:a,l 



rid) + tiand)cl da- + h.c] 
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al 
al 

- ^2[V la {2 - n d ) + Vi a {n d - l)]^^^,^, 



aa'l 



The first four te rms in ( |2.3| ) are the generalized Fano- 
Anderson model ( |2.2|) added by the Hubbard interaction, 
which is back to the Anderson model if t\ a = t\ a = t\. 
In the meanwhile, Via term takes into account the finite 
range interactions between the impurity orbital and con- 
duction electrons (I — 0,...,Nf). Without losing the gen- 
erality, we consider the screening channels as the spinlcss 
fcrmions and take Vi a = Via = Vi for I — l,..,Nf for 
convenience. The V X ai and V y i terms describe the spin 
exchange interactions. We take those with 1 = 1, —,Nf 
vanishing for convenience. The t2 term which is newly 
introduced in the Hamiltonian represents that two hy- 
bridizing channel electrons with opposite spins hop onto 
the impurity orbital simultaneously. To renormalize the 
model it is proper that this term as well as spin exchange 
terms are included in the Hamiltonian. 

In ( |2.3| ), if the coupling constants with the 'tilde' are 
equal to those without the 'tilde', the Hamiltonian is re- 
duced to the PH symmetric one under e d = —U/2. As 
we have shown that for this reduced model the consis- 
tency of the renormalization group equations forces the 
PH symmetry have to hold |23|] . This Hamiltonian also 
includes the Kondo model as a special case |23|. Spin- 
less models can also be reached by a proper choice of the 
coupling constants. 

We shall show there is a FL-NFL phase transition for 
the system. To understand the physical mechanics of 
the appearan ce o f the phase transition, we separate the 
Hamiltonian (2.3) into two parts, as done by Zhang, Yu 
and Su JT^] for the spinless impurity model, H — Hh+H s , 
where 

H h = ^2 e kc\ a c kl7 + t d n d + Un d ^n d i (2.4) 

k>0,a 

+ —r= y^[(*ig(l - nd) + l X an d )c\d a + h.c] 

V L a 

+ — y^X c i c -gdad-a + h.C.) 
a 

+ ^^2[V X<7 {2 - n d ) + V xa {n d - l)]4c a dld a 
+ T V y clc- a d f _ a d a 

+ j^2[V a{2 - n d ) + V a{n d - l)]4c CT 4' £ k'> 



J2 ekc li cki 

k,i>a 



1 x t 

+ — 2 yic\cin d . 

l>0 



(2.5) 



Here we have suppressed most of the suffices '0'. It is 
easy to see that the screening part, H s , of H is the 
multi-channel X-ray edge model, whose properties are 
controlled by a NFL fixed point; while the hybridizing 
part, Hh, of H is a one-channel impurity model, whose 
properties are given by some FL fixed points with dif- 
ferent symmetries. Since these two parts are coupled by 
the interactions between the conduction band and the im- 
purity orbital, there appears an interesting competition 
phenomenon between FL and NFL metallic behaviors, 
which causes the phase transition. 



III. PARTITION FUNCTION 

In this section we show that the partition function of 
the model at low energies may be equivalent to that of 
the multi-component kink plasma model with a weak ex- 
ternal field. 



A. Bosonization and Sum over Impurity Histories 

Since we deal with the low density system, the spec- 
trum of the conduction band can be taken as the linear 
form tk = (k — kf)/p with p = (hvp)^ 1 being the density 
of state at the Fermi surface. For the linear spectrum 
the conduction band can be treated by the bosonization 
technique || ^Gj. We first bosonize the screening part 
of the model, H s . We introduce the density operators: 

j W-k ^ W 

Pi( k ) = —ff c li c Q+k,i . Pi(- k ) = -7fJ2 c li c i- k > l > 

V L q=0 V h „-t 



q—k 



where W = X/tq is the conduction band width which is 
much greater than U and e d . The operators obey the 
Bose commutation relations 

\pi{k),pv{-k')] = kS k:k /Si tV , 

accordance with which, we can define standard bosonic 
creation and annihilation operators 

hi = -i\J^pi{k) , 4, = i^-pi(-k). 

The screening electron operator is represented, via the 
bosonic operator, as 



c\{x) = -=exp[-i<fo(x)], 



(3.1) 



and the screening Hamiltonian, H s , is transformed to 
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n. = E -olht - E 7=^ 

fej>0 ^ fc,Z>0 * 



L - bfeO^di ( 3 - 2 ) 



which can be diagonalized by the following canonical 
transformation 



T = cxp{i ^2 



pVi 



k,l>0 



{Ki + hi)n d }- 



Since this canonical transformation does not commute 
with d ai the hybridization terms are transformed to 
the coupling forms among the impurity, hybridizing and 
screening channels while other terms do not change in 
Hh- We have the effective Hamiltonian, H — Hq + Hj, 
with 



H o = E e k c la c k<y + Cd^d + U n d ^n di 

k>0,a 

+ j E^ 2 ~ nd ^ + Vv°( n d - l)Hc a di,d a , 

a. a' 



E -Am 



(3.3) 



k,l>0 ■ 



H, 



7= y^[(*i<r(l - n <i) + hand)^c\d a + h.c] 



J2(^ 2 cUl*d a d- a + h.c.) 



+ -fJ2°° C -' Td -° d < 7 - 
a 

The operator A* is given by A^ = exp{^ fc ;>0 -^!=L(i>£j- 



kL 

bki)} or in the coordinate representation A^(x) = 
exp{i^ — 4>i(x)} where Si = irpVi. 

TT 

Now, we tu rn to derive the partition function of the 
Hamiltonian (3.3), where H and Hj are considered as 
the unperturbed Hamiltonian and the interaction Hamil- 
tonian, respectively. Hence, we have a partition function 

Z = Trie- 13 ") = Tr(e- pH °T T exp{- / iJ 7 (r)dr}), 

Jo 

(3.4) 

where (3 is the inverse of the temperature, T T is the (imag- 
inary) time ordering symbol and Hj(t) = e H ° T Hie~ H ° T 
is the interaction-picture operator of Hi. Setting a se- 
quence of times < t\ < ... < r n < (3, the partition 

oo 

function is expanded as Z = Z n , where 



Z n = (-1)" [ dr n ... f 2 dT 1 Tr(e- f3Ho H I (T n )...H 1 
Jo Jo 



(n)). 
(3.5) 



The Hilbert space of H can be projected to four sub- 
spaces characterized by four impurity electron states, 
| a >= |0 >, \a > and |3 >= | jl>. Each term in the 
interaction H] plays a role of the 'dipole' operator caus- 
ing a transition between those sub-spaces. Such a kind of 
problems can be treated in Haldane's familiar way [j30| : 
Write the partition function in terms of a sum over histo- 
ries of the impurity. We will follow the formalism given 
by Si and Kotliar (T|j2|]. Setting P a = \a >< a\ is 
a projection operator to a-state subspace, Hq, then, is 
projected as 



Vi/ P 

7j 1± ct ± a.i 



where 



(3.6) 



(3.7) 



k>Q;a k,l>0 ^ 



The impurity levels, here, are projected as e d fl = 0, 
£d,cr = £d and £d,3 = 2e d ; the Hubbard interactions are 
U = U a = and U 3 = U. Making a shift of the 
ground state energy, we can redefine t d . a + U a as e a with 
£o = —(-d — U/4, e a = —U /4 and e 3 = e d + 3U /4, which 
obey ^2s a = 0. The potentials are given by 



Vr 



o -0, VJ, 



L 



V xa 



2V 0a , V 2 



xa 

17' 



Now, we sum over the impurity states by inserting a 
comple te s et of the impurity states at every discrete time. 
Thus, (3.5) is rewritten as 



i) n E / dTn - / dr 

n , „„ J0 Jo 



(3.8) 



ai } ai2 
^ 9a n -|- 1 a ?1 " ■ a 1 

xtr c {e- H ^+^-^0'(a n+1 ,a n )e- H ^ 
...e- H ^- T ^0'(a 2 , ai )e- H "i}, 



where tr c indicates the trace carried out for the conduc- 
tion band. a\, ...,ot n and n, ...,r n label a Feynman tra- 
jectory with Cin+i — Oil. g'a/3 are coupling constants, 

g'oAg'-vs) = 9iM<r) = hjVL(t la /VL), g '^_ a = g ' y = 

Vy/L and g' 03 = g' 2 = t 2 /L. The 'dipole' operators 
0'(a,j3) are given by 



O'(0,a) = O'\a,0) 
= 0'(-a,3) = 0' t (3,- ( r) = At4, 
O'(0,3) = O' t (3,0) = At 2 c{c|, 
0'(a, -a) = J_ a c a . 



(3.9) 
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A similar bosonization procedure to the screening 
channels can be applied to the impurity channel. As a 
result, the projected Hamiltonian transforms into 



H a =H c + e a - £ ^Jj(bl-b k .) 



(3.10) 



fc>0;p 



fc,Z>0 



where bka, the bosonic operator corresponding to Cka, is 
given through 

cl(x) = cxp[-i0 (T (x)], 



(3.11) 



The phase shifts 6° — irpV^L are defined by 

8° Q = 0, 8%, =6 + 8 xa 8 a , a , 8% = 2~S + ~S X<J , (3.12) 

where 5 _= npVoa, So = npVoa, <W = npV xa - and 
Sxtr = TtpVxo are the Born approximation values of the 
general phase shift s 2ta,rT 1 (iipVo ( j /2) , and so on. The 
Hamiltonian ( 3.1C| ) can also be diagonalized to H c + e a 
by the following canonical transformations 



U a = U a U s « = n CT exp{i<5^ CT }. 



(3.13) 



Meanwhile, the bosonized 'dipole' operators can be trans- 
formed to 



0(a,/3) = exp{iq^p(j) a }, 



(3.14) 



where <p a — {4> a ,4>i). The quantities which will be 
called 'charges', are 



( 1)<W H ) — , 

7T - - ' 

■$xcri 



(3.15) 



= (— -1)(<^_ CT -<W),0 

7T 



5o 



-(4a' 
7T 



28q - 26o), — 

7T 



%3 



5 xa > + 2(5p _ 25; 
7T 7T 



where a = (cr',l) and the 'charges' are anti-symmetric: 
= — gS. Using the bosonic 'dipole' operators 
( |3.14 ) and defining the interaction-picture operators 
0(a,/3)(r) = e HcT 0(a, (3)e~ H " T , the partition function 
reads 



(-ir E / ^ 







ro 



^ (3.16) 

T 



Tj+1 



x tr c (e~ ffc/3 0(a„ + i, a„)(T„)...0(a 2 , ai)(ri)), 



>Ei exp{- y^^+i 



r 



where g a p = g' a/3 To are dimensionless coupling constants 
and E a = e a To are regard ed as dimensionless external 
fields. The trace part of (3.16) is just the correlation 
function of n-vertex operators. In the zero temperature 
limit, it can be calculated by using the free boson field 
theory. The result is given rise to: 



/ / *H (3.17) 



X 



exp(-5[ri, ...,r n ;ai, ...,«„]), 



S[r,a] = E E C«, + i<-« J+1 ln " 

i<Cj a— (7,1 
- E ln (5a s + ia s ) + E ^i+l 



TO 

TO 



(3.18) 



which may be regarded to describe a one dimensional 
multi-component plasma of kinks with 'charges' q^g and 
'fugacities' g a p. E a is an 'external field'. 



B. Related to Charge-Spin Chain and Transitivity 
of 'Dipole' operators 



The partition function ( 3.17 ) describes a multi- 
component kink plasma with logarithmic interactions. 
There are six types of kinks which correspond to the 
different pairs of the impurity states, the charge pairs 
(Oct), (cr3) and (03) and the spin pair (f |). This leads to 
the logarithmic interactions between the kink-pairs can 
be transformed to the interactions among the charge- and 
spin-pairs. To reach this point we note that the 'charges' 
of kinks are of their transitivity: 

92/3 + ^ 7 = <&r (3-19) 
This leads to the rewriting of the action S[t, a], i. e. 

(3.20) 

TO 



S[t, a] = K(ai,a i+ i;acj,aj+i) ln ^~ 



i<j 



where 



K(ai,ai +1 ;aj,a j+ i) = K(a u aj) 



(3.21) 



+K(a i+ i, oij+i) - K(ai, a J+ i) - K(a i+ i,a>j), 

where K(a,(3), called charge- and spin- stiffness con- 
stants, are defined by [pl[ 



K(a,(3) 



1 

2 ^ 



(3.22) 



5 



The model described by ( 3.20| ) is equivalent to the one- 
dimensional charge-spin chain model with 1/t 2 interac- 
tion H] and incorporating the effective fields pf| : 

S[a,r] = ^2K(a l ,a J )- — -2— - ^ ln( 5cti+lQi ) (3.23) 



i<3 



Tj+l — Tj 
TO 



The model is just a special case of the general one- 
dimensional model with 1/r 2 interaction considered by 
Cardy p4| . It is worth to point out that the transitivity 
( |3.19 ) plays an important role to reach Cardy 's model. 
In fact, the transitivity of the 'charges' stems from the 
transitivity of the 'dipole' operators: 



0(a,p)0((3, 7 ) = 0(a, 7 ). 



(3.24) 



This transitivity implies that all 'dipole' operators to- 
gether with the vanishing operator <fi and the identity 
operator / form a semigroup under the operation: Ei- 
ther 0(a, /3)0(7, S) = <j) for /3 ^ 7 or the transitivity 
(3.24). Physically, the semigroup relates to the transi- 
tivity of the charge-spin- flips: e.g., the pair (03) can be 
separated as two pairs (Oct) and (ct3), etc. If the transi- 
tivity of the charge-spin-flips is broken, the charge-spin 
chain is interrupted. Therefore, the model is not self- 
consistent. Originally, the transitivity comes from the 
two-body hybridization a nd the spin-exchange interac- 
tion in the Hamiltonian ( |2.3j ). Therefore, the inclusion 
of them in (2^3) is proper to the self-consistency of the 
model. In the next section, we shall see that the cor- 
rect renormalization flows mustn't exclude those of the 
fugacities g-\i and goz- 

To conclude this section, we would like to mention that 
the transitivity of the 'charges' and their anti-symmetry 
leads to Jacobi identities, which implies the existence of 
a Lie algebra. We arrange an appendix ( appendix A) to 
discuss the properties of the Lie algebra. 



IV. PHASE STRUCTURES 

In this section, we study the phase structures of the 
models by using the AYHC poor man's scaling theory. 
The details of the derivation of the renormalization group 
equations have been given in Appendix B. Using those 
equations, we analyze the phase structures of the general 
model and various reduced models. 



A. General Model 

In the most general case, six non-zero fugacities g a p 
are independent of each other. Furthermore, six stiffness 
constants are the functions of six phase shifts Sq, Sq, 8 xa 
and 8 xa in the quadrate. The six stiffness constants are 



also independent of each other. To see thi s poi nt, we con- 
sider the small phase shifts. In this case ( ft.22| ) expresses 
a non-singular linear transformation if the quadrate are 
neglected. Thus, the renormalization group equations 
are exactly the 4-state example of the Cardy's general 
results |24[| and their symmetry breaking generalization 
p5[ (also see Appendix B): 



dgi 



dlnTo 

dg v 

dlnTo 

dgia 
dlnTo 

dgi 
dlnTo 



Ko a )g llJ + g l ,- a g y e- E -^ E « +E ^ 

E 3 + \(E +E a ) 



(4.1) 



#23l-<re 

(1 - K a - a )g y + giagi-ae 
h<rh,-*e- E ° + H E ° +E -°\, 

(1 - K-a3)gia + 9\-cr9lZ 
,-E, + \(E 3 +E- a ) 



-Eo + ^Ez+E-v) 



9 v 9i, 

(1 - ^03)32 + } ,fllr.7l- - < 



for the fugacities. Here K a p = —K(a,{3). For the exter- 
nal fields, i. e. U and e<j, one has 



to 



de d 
dlnTo 



rotd + (<7i T + 9 2 n )e Eo 
1 



En —E"t 



(4.2) 



1 



n 2 \ P E r~ E 3 



To 



dU 
dlnTo 



t U + ( 5 2 T + g 2 n ){e l 



9le E ° 



- E3 



.Eq — E-\- 



A —E 3 



- gl{e Ea ~ E3 + e E ^ Ea ). 
For the stiffness constants, one has 

z E f>- E i{K af) + K^-K ai ) 



dK af} 
dlnTo 



K 



Pi) 



(4.3) 



7 

7 



We make a few comments on the renormalization group 



equations (4.1), (4.2) and (4.3) 



(i) As we have seen, if the model can be mapped to the 
special model of Car dy's, V y and ti can not vanish. Fur- 
thermore, from (4.1), we see that even we start from the 
vanishing initial values of gi and g y , they will increase as 
g\ a and g\„. This means that the flows of the vanishing 
values of them are not correct renormalization flows. 

(ii) Since our renormalization group is perturbative 
in its treatment of g a p, the phase structure of the sys- 
tem can be described by g a p. The strong-coupling fixed 
points in the quantum impurity problem are outside the 
range of validity of the AYHC poor man's scaling theory, 
but it can give the correct flow directions. In the strongly 



G 



coupling limit, the fixed point behaviors need to be more 
carefully discussed. 

(iii) Near the fixed point, e d , U — > in the order of 
9a/3- Therefore, we can regard all e Ea = 1 when we dis- 
cuss the behaviors near the fixed point. Then, we can 
rewrite the renormalization group equations in the form 
of perturbative g a p as follows: 



dgia 
dlnTQ 

dg v 

dlnTQ 

dlnTQ 
dg2 
dlnTQ 



(1 - K 0(T )g la - + gi.- a g y + gsh-v, ( 4 -4) 
(1 - K a - a )g y + g\ a g\- a + gi a gi-a, 
(1 - K^vzjgxv + gi-crg 2 + g y g\-a, 
(1 - K Q3 )g 2 + ^2 9i<t9i-<t, 



TO 



TO 



(fed 

dlnTQ 

dU 
dlnTQ 



= T ed - 



I 2 



9 2 n 



9\\ 



(4.5) 



~2 \ i 2 2 

ffll) +52 ~9y 



TQU + 2(gl-gl) 
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Under those considerations and according to (4.4), we 
can describe the phase structure of the system. 

(i) For all K a p > 1, all fugacities are irrelevant and 
renormalized to zero. There exists a weak coupling fixed 
point <7*^ = 0. The fixed point Hamiltonian is similar 
to that of the multi-channel X-ray edge problem [p2[ . 
The system exhibits a power law decay of the correla- 
tion function with a non- universal exponent. This is a 
NFL phase. By the definitions of Ko a , the existence of 
the screening channels ensures the possibility of K 0a > 1 
and then reaches the NFL phase. The one-channel model 
in a small exchange interaction can not have this phase 

(ii) Once one of the coupling constant becomes rele- 
vant, the low temperature fixed point of the system, in 
general, is Fermi-liquid-like. It is still possible that there 
is a plenty of phase structures at finite temperature. Be- 
cause all K a p are independent of each other, in general 
there will be a unique smallest K a p with (a/3) = (ao/3o)- 
When K ao p decreases to 1 and then smaller than 1, a 
phase transition happens. There are such six phases that 
5q /3o i s relevant and all the g a $ ^ 3a /3 are irrelevant in 
each phase. They may have different finite temperature 
behaviors. For example, we consider (^-relevant phase. 
We focus on the fixed point Hamiltonian. Introducing 
the scaled parameters V y /L = Jj_/2, Vq /L = —J z /4: 



and V xa /L — J z /2 with J z > (this ensures < 1), 
we see the fixed point Hamiltonian H* is just the one- 
channel Kondo Hamiltonian: 



H* = ^2 e kcl„Ckcr + e d n d + Un d ^n d i 



(4.7) 



k>0,a 



+ E -Ab kl + ^{sts^+h.c.)+j z s z c si 

k>o,l = l P 

Here S c and s,j are the spin operators of the hybridizing 
and local electrons respectively. The g y ( JjJ-relevant 
phase is controlled by the strong fixed point of Kondo 
problem. Furthermore, the values of Vq g and V xa are 
not important as long as they are within the parameter 
regime of the ^-relevant phase. Hence, the fixed point 
corresponding to (^-relevant phase is the FL, which is 
a collective mode of the system and transfers to a state 
with a local moment at T > Tk, the Kondo temperature. 
Equations (4.5) is reduced to 



TQ 



de d 
dlnTQ 



T e d 



■9y, 



TQ 



dU 
dlnTQ 



TqU + 2g 



which imply that U increase positively if its initial value 
is not negative. This repeats the results of Anderson 
model: When the impurity level width T *~ t\ is much 
less than U, the model is equivalent to the Kondo model 

For the gi cr (gi (T )-relevant phase, according to (4.5), U 



keeps to be almost invariant if its initial value is the same 
as at the fixed point. Hence, r 3> U in this case and the 
groundstate of the system may be read off from the ex- 
actly solvable Fano- Anderson model. It is a FL but there 
is no a local moment in excitation states. So we call this 
phase the free-orbital phase Q . 

For the 52-relevant phase, the empty and doubly- 
occupied states favor over singly-occupied states. The 
fixed point Hamiltonian H* = Hq + Hj where Hq can 
be read off from (3.3) with replacing the all of param- 
eters by their fixed point values and Hj has a leading 
term R ~ (g^/To) ^2 a -{ c t c -<jdad- (7 + h.c). The matrix 
elements of R in the fixed point can be evaluated by a 
similar considerations in the X-ray edge problem fl33| . 
As done in |Tq|, define \n d — > to be the eigenstates of 



i?o|n d =o and \n d = 2 > to be H*\ nd =2S. Then 
< n d = Q\R\n d = 2 >~ r Q . 



(4.8) 



The anomalous exponent a = (1 — Kqs)/2. (If a < the 
operator R is irrelevant, and the fixed point is the NFL 
as we have discussed.) In the present parameter regime, 
a > and the operator R are relevant. The fixed point 
is regarded as the FL |T^| . 

Thus, all phases with one relevant fugacity are FL but 
they may have different finite temperature behaviors. A 
FL-NFL phase transition occurs as the model parameters 
vary. 

(iii) When several of the K a p < 1, the phase transi- 
tion can be characterized by Cardy's graphical method 
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J24j . Draw a graph whose vertices are labeled by the 
states 0,cr and 3 (see Figure 1(a)). A particular tran- 
sition corresponds to drawing a set of the edges (a/3) 
which express that those K a p < 1. When K a p and Kpi 7 
become smaller than one, g a p and gpi-y are relevant. If 
f3 7^ /?', g ai is still irrelevant as long as K ai > 1. Before 
the next phase transition happens, we change the edge 
linking the vertices a and (3 to a new vertex and the edge 
of (3' and 7 to another vertex. Then repeat the previous 
process when other K a > 7 > < 1. If/3 = (3', even K ai > 1 
the second and third terms of ( }4.1| ) also drive g aj to be 
relevant. Reflecting this point graphically, the vertices 
a, /3 and 7 are linked as an edge which is changed to 
a new vertex before the next phase transition appears. 
Cardy had explained the physical reason to cause g a /3 
becomes relevant p4| : Logarithmic interactions between 
kinks are no longer existent in the linked piece of the 
graph. Repeat the process till all the vertices are linked 
so that no further transitions are possible. Cardy calls 
the final phase high temperature (HT) phase, in which all 
the fugacities are relevant. Notice that the renormaliza- 
tion group equations are not valid quantitatively in the 
strongly coupling region. Therefore, the phase structure 
of the model in the relevant phase must be discussed more 
carefully. In some situations, the real phase diagram of 
the relevant phases is different from that given by the 
graphical method. To see this we now turn to various 
reduced models from the general model. 



B. Spin Symmetric Model 

The spin symmetric model is a model with g a ^ = g a i 
and K a ^ — K a ^. The phase structure of the system 
can be described graphically in Figure 1. Every diagram 
row of Figure 1 represents a sequence of transitions from 
the NFL phase to HT phase. The diagram consisting 
only of the vertices represents a phase while the diagram 
with lines is a critical point. For example, Figure 1(a) 
represents the sequence of transitions corresponding to 
K a < K 03 ,K 3a ,K n . In fact, Fig. 1(a) and 1(b), 1(c) 
and 1(d) reflect similar physical processes, respectively. 
The four-vertex phase in the figure is the NFL phase. 
The two-vertex phase of l(a)(l(b)) is the girr(gi,<j) and 
g v relevant phase. We arrive at the strong-coupling fixed 
point of the asymmetric Anderson model. This is the 
FL. The three- vertex phase of 1(c) (1(d)) is controlled 
by the strong-coupling fixed point of the Kondo problem. 
And two- vertex phase is the same phase as 1(a) (1(b)). 
The three- vertex phase of 1(e) is the (72-relevant phase, 
which is the FL. The phase structure of the spin sym- 
metric model is different from that of the general model: 
The free orbital phase does not exist. Physically, since 
9a] = 9ai the spin can not be frozen out. Furthermore, 
when the (0 \) kinks are bounded, the (0 J.) kinks are 
also bounded because of the spin symmetry, which leads 
to the (tl) kinks are bounded too. Therefore, g y be- 



comes relevant together with g\ a or g\ a . The free orbital 
condition, r > P, can not be satisfied in the spin sym- 
metric model. We would like to mention a point that the 
spin symmetry does not change the dimensions of the Lie 
algebra described in Appendix A. Correspondingly, the 
largest number of the critical points in the sequences of 
the phase transition keeps to be invariant. As we have 
mentioned, in the strong coupling limit, the graphical 
method is not always valid quantitatively. Here this point 
reflected in the phase diagram is that the phase transi- 
tions from three-vertex phases to two-vertex phases do 
not really happen because both of the asymmetric Ander- 
son model and the Kondo model have the same strongly 
coupling fixed point. 



C. Particle-Hole Symmetric Models 

Setting g\ a = <7i j(T and Ko a — K- a3 in the general 
model, the consistency of the renormalization flows of g\ a 
and gi a forces that E 3 = E (as long as their initial val- 
ues are equal). Thus, if the PH symmetry of the common 
hybridization term in the model is not broken, the con- 
sistency of the renormalization group equations requires 
that the system keeps the symmetry. The phase struc- 
ture of the PH symmetric model is described in Figure 
2. Since there is no spin symmetry, the free orbital phase 
exists (the two-vertex phase in Figs. 2(a) and 2(b)). The 
phase described by two-vertex diagrams in Figure 1 dis- 
appears because the model is PH symmetric. The Kondo 
phase and <?2-relevant phase still exist. 

If we consider spin and PH symmetric model, the phase 
structure is shown by Figures 3(a),(b) and (c). This 
model has been discussed in Q. Figures 3(a), (b) and 
(c) can be reexpressed more explicitly by Figure 3(d). 
Notice that there is no three critical point sequence of the 
transition in the PH symmetric model because it is im- 
possible that K^i < Kq 3 < K 0a and K 03 < K n < K 0a . 
Physically, those sequences of transition break the PH 
symmetry and then are forbidden. Relating to the Lie 
algebra of the plasma 'charges', under the PH symmetry, 
the dimensions of the Lie algebra reduce to 3. As we 
have done the strongly coupling limit needs to be 
treated more carefully. The real phase diagram is shown 
in Figure 3(e). The phase of all fugacity relevance is di- 
vided into three parts: / is still in the Kondo strongly 
coupling phase; 77 is the same as GVphase in which the 
empty and doubly occupied states are favored; and III 
is all impurity state mixed phase. 



V. MIXED VALENCE STATE REGIME AND 
FIXED POINTS 

In the previous section, we term the FL and NFL 
states only in the sense that whether the infrared or- 
thogonality catastrophe of the X-ray edge type appears. 
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In fact, the situation may be subtler if the we consider 
the local state properties in the problem Jl^,^|. To see 
this, we study the fixed points in the mixed valence state 
regimes, in which the renormalizcd impurity level shift e* d 
is of the same order as the renormalized resonant width 
T* = 7T5* 2 /t, i. e., \e* d \ < V*. Most interesting phys- 
ical phenomena appear in the particle-hole symmetric 
model. At finite U, the particle-hole symmetry implies 
e c i = —U/2 for their bare values. The renormalization 
group equations show that this constraint is kept after 
the renormalization except in the case of infinite U. So, 
the crossover between various regimes is still determined 
by the ratio |e^|/r*. We can still define the mixed valence 
states if |e^| ~ T* . On the other hand, the particle-hole 
symmetry reduces the number of the independent stiff- 
ness constants to two , 70 and 7^. The fixed points in the 
7o — "fx plane can be figured out according to the stiffness 
constants e tl = ±70, e x = -f x e t2 = 270 - J x - The weak 
coupling mixed valence states is within the renormalized 
parameter regime el > 1 , e* > 1 and e| 2 > 1. The 
strong coupling mixed valence states are corresponding 
to < 1 and e£ < 1. There are two intermediate cou- 
pling mixed valence state regimes in which the renormal- 
ized stiffness constants are given by 1) e£ < 1, e* < 1 
and e* 2 > 1 and 2) e* ti > 1, e* < 1 and e t * 2 < 1. The 
latter associates to the intermediate coupling mixed va- 
lence state we have already known in the infinite U case 
po| , ^5| while the former to the G2-irrelevant mixed va- 
lence states we newly find in the finite U case here. We 
will not repeat the fixed point behaviors in the strong 
coupling case and the known intermediate coupling case 
but refer to the interested readers to read the relevant 
literatures |3(],^5| . What we would like to do is to review 
the weak coupling states and focus on the G2-irrelevant 
states. 

Within the weak coupling mixed valence regime, all 
running fugacities could be obtained by solving the lin- 
ear renormalization group equations. As we have seen 
in the last section, there is a NFL fixed point in the 
sense of the X-ray edge infrared orthogonality catastro- 
phe. Furthermore, in the mixed valence regime, there 
are also the violation of the Fermi liquid behavior in the 
calculation of the local single-particle and two-particle 
correlation functions. Some of them have been shown in 
refs. jl^m such as Grf c , Gdd and two-particle excitonic 
and d-electron correlation functions as well as dc correla- 
tion functions. Divergences appear in the Gd c , Gdd and 
the dc correlation functions, which shows the local NFL 
behavior. At the finite U , one may have more divergent 
correlation functions. For instance, the single d-electron 
paring function 

Gd-pair(r) =< 0K(rK,(0)|3 >~ r~ a \ 
or 

G d _ pair (o;) ~ u~ 1+a \ (5.1) 
with the exponent 

a d = 26* 2 -5* x 2 (l-S a ^). 



The two-particle d-electron paring correlation function is 
given rise to 



<d a d a ,( T )didUo) >~r- a \ 

with the exponent 



(5.2) 



CC4 = 



6* + « 



■2\ 2 



We see that both d-electron pairing correlation func- 
tions are divergent in the low energy limit. 

The G2-irrelevant intermediate coupling mixed valence 
state is newly discovered at the finite U model. The 
renormalization group equations yield an initially de- 
creasing the paring charge fugacity 172 and initially in- 
creasing the charge and spin fugacities g\ and g x . It can 
occur for which the spin and charge kinks are unbound 
while the paring charge kinks are bound. The two-body 
hybridization is irrelevant while the hybridization and 
exchange coupling are relevant. Therefor e, th e correla- 
tion functions such as < 0|i?|3 > given by (4.8) and local 
pairing single- and two-d-electron correlation functions 
are still divergent. This gives a new set of local NFL 
fixed points. 



VI. CONCLUSIONS 



We have constructed a general impurity model with 
two-body hybridization. Through comparing with the 
partition functions, via bosonization, we have found that 
this impurity problem is asymptotically equivalent to 
a six-component plasma of kinks with external fields. 
Hence, the renormalization group analysis can be sys- 
tematically carried out in the framework of the AYHC 
poor man's scaling theory. We emphasize that it is im- 
portant to introduce the two-body hybridization in the 
Hamiltonian. This ensures the transitivity of the 'dipole' 
operators and then the equivalence of the partition func- 
tions. Moreover, the renormalization group equations 
show that the flows excluding the two-body hybridization 
are not correct renormalization flows. In some phase of 
the phase diagram, 172 is relevant, i.e., hybridizing elec- 
trons hop pairly onto and off the impurity orbital . This 
behavior of the electrons resembles the 'Cooper pair' in 
BCS theory. The property of the (^-relevant phase is 
not clear yet. Perhaps, it relates to the superconduct 
phase in an impurity lattice model. In the phase dia- 
gram, there is a NFL phase which controlled by the fixed 
point behavior of the multi-channel X-ray edge problem. 
Around this phase we have six other phases called the 
(72-relevant phase, the strong-coupling Kondo phase and 
four equivalent free orbital metallic phases. Hence, there 
is a FL-NFL phase transition in the system. There are 
six sequences of phase transitions from the NFL phase to 







the HT phase in which all fugacities are relevant. Each 
sequence has three critical points and four phases, which 
is same as the dimensions of the Lie algebra correspond- 
ing to the plasma 'charges'. 

We have discussed various reduced models. For the 
spin symmetric model, the free orbital states do not exist 
because we can not freeze only one spin state in the impu- 
rity orbital and the free orbital condition can not be sat- 
isfied. For the PH symmetric model, we addressed that 
the PH symmetry of the common one-body hybridization 
term in the Hamiltonian forces = — U/2 by the con- 
sistency of the renormalization group equations and also 
decrease a critical point in the sequence of transitions, 
correspondingly reducing a dimension of the 'charges' Lie 
algebra. 

We also discussed the fixed point behaviors in the 
mixed valence state regime. Except the X-ray edge type 
infrared orthogonality catastrophe in the screening chan- 
nel, we also found the local NFL signal in the weak cou- 
pling and intermediate coupling mixed valence states. 
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APPENDIX A: LIE ALGEBRA AND ITS 
RELATION TO PHASE STRUCTURES 



(ii) The Lie algebra is 4-dimensional, which is same 
as the dimensions of the impurity Hilbert space. The 
structure constants form the adjoint representation by 
(T a )p~ ( — cp a a. For the general parameters 5q, 5q, S xa and 

8 xa , aTo + bT-\ + cT± + dTj, = only if all coefficients a, 6, c 
and d are zero. 

(iii) In the infinite U limit, the doubly-occupied state 
is completely suppressed in the discussion of the low 
energy behavior of the system. The physics has been 
discussed in |Q [ [LslJ^ . The impurity Hilbert space is 
3-dimensional, and correspondingly, the Lie algebra re- 
duces to 3-dimensional one: £ — {To, T a }. 

(iv) In the PH symmetric model, q^ a — q~L a3 , which 
leads to the dimensions of the Lie algebra is reduced to 
3 too. This consequence can be checked in the adjoint 
representation by the following linear transformation 

L =T + T 3 + T t + T u (A3) 
L 1 =T -n + T^+T i , 
L2 = T — T 3 + Tf — Ti , 
L3 = T + T 3 — Tj — Ti . 

The transformation is a non-singular but one finds that 
L3 = 0. This implies that the dimensions of the Lie 
algebra is 3. 

At the moment, it is not clear that what are the phys- 
ical implication of the Lie algebra and the corresponding 
Lie group. We only know that there is a correspondence 
between the phase structures and the dimensions, D, of 
the Lie algebra: D — 1 is equal to the largest numbers 
of the critical lines between the NFL phase and the HT 
phase. 



In this appendix, we discuss some mathematical prop- 
erties of the 'charges' g"^. The screening components, 
q a g, are trivial and we focus on the spin-components, 
<f a p. We define q^p are given by the spin-components for 
7 = a or va nish for 7 = and 3. In terms of the transi- 
tivity ( 3.19 ) and the anti-symmetry of the 'charges', we 
immediately have q 1 ^ obey the Jacobi identities: 



APPENDIX B: 



lipids 

7 



0. 



(Al) 



which means that the 'charges' may be regarded as a set 
of structure constants of some Lie algebra. Denoting the 
generators of the Lie algebra by T a , then wc have 

[To, T a ] = ql a T a , , [Ta, T a ,} = qCT a „ , (A2) 
[T CT ,T 3 ] = q^' 3 T a >, [T ,T 3 ] = q^T*. 

The Lie algebra has the following properties: 

(i) It is solvable. We call the Lie algebra £. Equa- 
tions (A2) imply that [£,£] C £ lt [£ 1 ,£ 1 ] C £ 2 and 
[£2, £2} = 0, where £1 = {T a } and £ 2 = {q^T a }. This 
means that £ is solvable. 



RENORMALIZATION GROUP 
EQUATIONS 



Our partition function (3.17) has arrived at a special 
form of the Cardy's general model with an external field 
added. Therefore, the renormalization group equations 
can be directly derived by using the AYHC poor man's 
scaling theory. The details of the derivation have been 
presented in |2^] 



[18 



|25|] . Here we just address some- 
thing in their derivation. 

The renormalization group equations describe the 
flows of the dimensionless couplings, the fugacities g a @, 
the stiffness constants K(a,(3), and the external fields 
E a , as the band width is reduced. 

Firstly, when the cut-off To — > tq + dro, to compensate 
the changes of to's power terms and external field terms 
in the partition function ( |3.17 ), the fugacities and the 



external fields change as follows: 

= (dln,T )g ai+iat (l + K(a i+ i,ai)), 
SE ai — E ai dlnTQ. 



(Bl) 



The cut-off of integral limits also causes the variation 
of the logarithmic terms. After integrating out the Tj+i, 
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we obtain the variation involving two such kinks i 5 i + 1 
and a third kink j which is given by 



exp{K(ai, a,; + i; otj, aj + i)ln— 



to 

T, - Ti 



TO 



TO 



+K(a i+ i,a i+2 ] aj,<Xj + i)ln 

exp{E ai+2 - E ai+1 } 

PS (dlriTo)g ai a i+1 9a i+1 a i+2 

exp{(K(ai, otj) + K(a i+ 2,a j+ i) 

-K(a h a j+1 ) - K(a i+ 2,aj))ln— — — } 



TO 



exp E aiJ ^ 2 E ai+1 . 



(B2) 



The only approximation is Tj — Ti To, which means the 
kinks are rare enough. 

As long as on 7^ a>i+2, this can be incorporated into a 
renormalization of g ai a i+1 , i.e- 



(B3) 



after integrated over the 'closed pair' Q range r,_i < 
Ti < Ti + 2 • This can be interpreted as the renormalization 
of the stiffness constants: 



dK{a,p) 
dlnro 



(A-(o,/3) + A-(a, 7 )- K((3,j)) 



(B8) 



Ea—E-y 



(K(a,0)+K{p,i)-K{a,i)). 



Equations (B4)-(B6) and (B8) are our renormalization 
group equations in the model. The relations g a ,a = 
and K(a, a) = and £" Q = are preserved in the 
renormalization process. 



which together with (B.l) and g a p = gp a leads to the 
renormalization flows of the fugacities in arbitrary exter- 
nal fields: 

Jg = (l + K^mg^ + ^g^pe-^^l 

1 

(B4) 

However, for the neutral pair with on = the lead- 

ing term of the left hand side of (B.2) is independent of 
the relative position of the kinks. It becomes: 

/ 7 j n 2 T i+2 — Ti-i E E 

(dlnT )g a e - -h-i. 

T 

which leads to the renormalization of the external fields: 



5E n 



-dlnro^glpt 



E a -E a 



Taking into account ^2 E a — and the renormaliza- 
tion of the free energy F , the renormalization flows of 
E a and F are 



dE a \ - 2 

7 

dFrp _ _ 1 \p 2 

ap 



E a —Eq 
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, (B5) 



(B6) 



Then, we must consider the next term in the expansion 
of (B.2) as to/(tj — Ti), which is 

(dlnTo)gl. a . +1 K(ai,ai +1 ;aj,aj + i)(ln(Tj - r. i+2 ) (B7) 
-ln(Tj - Ti_i))e Bcv ' _B °'+ 1 , 
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Figure Captions 

1. The phase structure of the spin symmetric 
model, (a) K 0a < K a3 ,K 03 ,K n . (b) K a3 < 
K aa ,K 03 ,K n . (c) K n < K 0a < K a3 ,K 03 .(d) 
K n < K a < K 0rT ,K 03 . (e) K 03 < K 0a ,K a3 ,K n . 

2. The phase structure of the PH symmetric model, 
(a) Kot < K 0l ,K n ,K 03 ; (b) K 0L < K 0h K n ,K 03 : 
{c)K n < K 0rT < K 03 ; (d) K 03 < K Qa < K n or 
^03 < K 0a < K n < K - a . 



Fig. 3. The phase structure of the spin and PH sym- 
metric model, (a) K 0a = = K n3 ; (b)Kn < 
K a < K 03 ; (c) K 03 <K 0<7 < K n . (d) The phase 
diagram in 70 — -f x space, where 70 = (1 — 5 x /ir — 

<y^ 2 + (<V^) 2 + Ez(<V^) 2 and 7* = (i-4A) 2 . 
(e) The real phase diagram in 70 — "f x space. The 
thick lines are phase boundaries. The thin lines 
divide the Fermi liquid phase into the different re- 
gions characterized by different behaviors at finite 
temperature. 
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